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1.  Introduction 


Because  dielectric  waveguides  are  frequently  used  as  part  of  the  ex¬ 
perimental  setups  in  the  laboratory,  it  is  of  interest  to  have  general  pro¬ 
grams  that  calculate  their  characteristics.  These  waveguides  generally 
consist  of  thin  sheets  (a  few  wavelengths  or  less)  of  a  low-loss  nonmag¬ 
netic  dielectric  material.  The  electromagnetic  wave  is  launched  at  one  end 
of  the  slab  and  guided  to  the  sample  [  1 ). 

In  this  report,  we  idealize  the  situation  to  a  single  slab  of  infinite  length  in 
two  dimensions,  and  a  thickness,  2d,  which  will  be  arbitrary.  First,  we  in¬ 
troduce  the  appropriate  Maxwell’s  equations  for  the  problem  and  derive 
the  wave  equation  for  the  transverse  electric  (TE)  and  transverse  magnetic 
(TM)  modes.  We  obtain  the  general  solutions  for  both  modes  inside  and 
outside  the  dielectric  and,  by  satisfying  the  boundary  conditions  at  the  sur¬ 
faces  of  the  slab,  obtain  four  transcendental  equations  which  provide  the 
solutions  to  the  problem.  There  is  one  equation  for  each  mode,  the  TE 
even,  TE  odd,  TM  even,  and  TM  odd.  Next,  we  assign  values  to  our 
parameters  and  numerically  solve  the  transcendental  equations.  This  pro¬ 
cedure  is  repeated  for  several  thicknesses  and  indices  of  refraction.  The 
results  are  presented  in  graphical  and  tabular  form.  The  appendix  contains 
the  listings  of  the  programs  used. 
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2.  Maxwell’s  Equations 


Our  waveguide  extends  infinitely  in  the  y  and  z  directions,  and  from  -d  to 
d  in  the  jc  direction  tsee  fig.  la).  The  wave  is  launched  in  the  y  direction. 
The  relevant  Maxwell  equations  for  the  problem  (in  centimeter-gram- 
seconds)  are 


V  x  £  =  ikH 


(1) 


and 


-V  x  H  =  ikeE  ,  (2) 

where  e  is  the  dielectric  constant.  The  time  dependence  is  assumed  to  be 
exp(-icot),  with  k  =  o )/c.  Throughout  the  problem,  (i  is  considered  con¬ 
stant  and  equal  to  one.  In  all  cases  considered,  we  assume  that  none  of  the 
fields  vary  in  the  z  direction  (d/dz  =  0). 

In  our  first  case,  the  TE  mode,  £  =  0,  £  =  0,  and  £,  *  0.  So,  Maxwell’s 
equations  yield 


ikH,  =  E,  , 

(3) 

(4) 

and 


a 


Hx 


=  ikzEz  . 


Combining  these  equations  gives  us  the  wave  equation 


V,2  Ez  =  -k2  eE2 


where  V?  =  ^ 

oxr  a.v 


e  =  e 

e  =  1  [<\>d  , 


(5) 


(6) 


« 


Figure  1.  Physical  repre¬ 
sentation  of  problem 
with  (a)  transverse 
electric  field  and 
(b)  transverse  magnetic 
field. 
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We  will  assume  a  dependence  in  the  y  direction  of  the  form  exp(t'p.y).  We 
will  also  assume  odd  and  even  solutions  of  the  form  sin  kxx  and  cos  kxx 
inside  the  dielectric,  and  a  decaying  exponential,  exp(-rr),  outside  the 
dielectric.  For  guided  waves  in  a  lossless  material,  we  must  have  T  real 
and  positive. 


For  even  solutions,  using  equation  (6),  we  have 
Ez  =  E\  cos  k\xel^y 


Ez  =E2e~rxe‘to 


\x\<d 

\x\>d 


where 


therefore. 


k]  +  [32  =  k?e 
-r2  +(32  =k2  ; 


r  =  V*2(e- \)-k] 


We  know  that  Ez  and  Hy  must  be  continuous  at  the  boundaries  x  =  d  and 
x  =  -d.  We  need  only  solve  at  one  of  these  boundaries,  though,  because 
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of  the  symmetry  of  the  problem.  We  can  find  Hy  by  using  equation  (4) 
with  equations  (7)  and  (8): 


Hy  =  E\  k i  sin  k\  xe‘^y 
Hy  =  ~  ve2 <rTxeh 


Ul<t/  , 


Setting  equation  (7)  equal  to  (8),  and  (12)  equal  to  (13)  gives  us  the  sys¬ 
tem  of  equations 


which  is  reduced  to 


E\  cos  k\  d  =  E2e'Yd  , 
E\k\  sin  k\  d  =  E2Te~Yd  , 


k\  l2Lf\k\d  =  V  . 


Throughout  the  remainder  of  our  discussion,  we  will  be  using  the  follow¬ 
ing  substitutions: 

i<  =  M  (17) 


A2  =  k2(e-\)d2  . 


Using  equations  (11),  (17),  and  (18),  equation  (16)  becomes 


tan  u  = 


This  transcendental  equation  will  be  solved  in  section  3. 

For  odd  solutions,  sine  replaces  cosine  in  the  material,  making  our  wave 
equation  reduce  to 


Ez  =  £3  sin  k\xe‘Py 


fvl  <  d  , 


m 


Equations  (9),  (10),  and  (11)  still  hold.  Using  equation  (4)  on  equations 
(20)  and  (21)  yields 


Hy  =-jj^E-jk\  cos  k\xel^y 


Hy  =jjcE4re~rxe‘^ 


Ul  <  d  , 

(22) 

bc|>d  . 

(23) 

At  x  =  d,  we  have  a  system  of  equations  which  produces 

T d 


-cot  k\d  - 


kid  ■ 


Using  our  substitutions,  from  equations  (11),  (17),  and  (18)  we  get 


-cot  u  = 


iA2^2 


(24) 


(25) 


This  transcendental  equation  will  also  be  solved  in  section  3. 


The  calculations  for  the  TM  modes  are  very  similar.  We  use  the  same 
Maxwell’s  equations,  but  now  the  electromagnetic  wave  is  described  by 
figure  lb: 


Hx  =  0,  Hy=Q,  Hz*0  . 


Maxwell’s  equations  yield 


F  =-—  d  H 

x  ike  <5y  z  ’ 

(26) 

F  =—  d  H 

y  ike  z  ’ 

(27) 

Hz=Tk  ^Ey~^Ex)  • 

(28) 

Combining  these  equation^  yields  the  wave  equation 

Vf  H z  =  -k2  eHz  . 

(29) 

For  even  solutions,  we  assume 

Ht  =  H\  cos  k\xel^y 

M  <  d  , 

(30) 

Ht  =H2  e'rxeifiy 

[x\>d  , 

(31) 
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where  and  T  are  given  in  equations  (9)  and  (10).  E  and  Hz  must  be 
continuous  over.t.  Using  equation  (27)  with  equations  (30)  and  (31)  gives 
us  a  system  of  equations  at  x  =  d  which  yields 

tan  k\d  =  .  (32) 

Or,  after  making  our  usual  substitutions,  with  equations  (11),  (17),  and 
(18),  we  get 


tan  u  =  e 


lA2  -  u 2 


This  will  be  numerically  solved  in  section  3  with  the  other  transcendental 
equations. 


For  the  TM  odd  mode,  a  very  similar  set  of  calculations  arrives  at  the 
solution 


-cot  ll  =  £ 


JA2-u2 


We  now  have  the  four  transcendental  equations  providing  the  even  and 
odd  solutions  to  the  TE  and  TM  modes. 
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3.  Numerical  Calculation  of  Solutions 


In  section  2  we  obtained  analytical  solutions  to  our  problem  in  the  form  of 
these  four  transcendental  equations: 


TE  even 


TE  odd 


TM  even 


TM  odd 


tan  u  = 


-cot  it  = 


JA2  -u2 


M2  -  lir 


tan  u  =  e 


-cot  u  =  e 


M2  -  u2 


M2  -  M2 


We  wish  to  numerically  calculate  all  the  solutions  to  these  equations  using 
the  Newton-Raphson  method.  We  must  first  assign  values  to  A  and  e.  To 
assign  a  value  to  A  we  must  first  break  it  down  to  its  components.  We  as¬ 
sume  the  following: 

n  -  3.5  (GaAs) 

8  =  12.25  =  n2  (GaAs) 
d  =  0.175  (am 
X  =  1 .064  (am 
k  =  5.93  x  104  (am 
A  =  3.466  =  -Jk2  (e  -  1  )d2  . 

These  values  are  chosen  to  illustrate  the  method  of  solution  to  be 
described,  but  are  not  to  represent  practical  values  in  use  in  the  laboratory 
|2|.  Although  there  may  be  some  use  of  waveguides  as  thin  as  this,  most 
optical  waveguides  are  in  the  neighborhood  of  1  (am. 

Just  knowing  these  parameters  is  not  enough.  We  must  know  the  number 
of  solutions,  and  their  approximate  values,  for  the  Newton-Raphson 
method  to  work  well.  Figure  2a  is  a  plot  of  tan  u,  -cot  u,  (A2  -  u2)m/u, 
and  e(A2  -  u2)]l2/u  versus  u.  The  intersections  of  the  trigonometric  func- 
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tions  with  the  latter  two  functions  are  the  solutions  to  our  Droblem.  Their 
exact  values,  so  far,  are  unknown,  but  they  provide  a  good  first  guess  for 
the  program,  as  well  as  a  check  on  the  number  and  range  of  our  solutions. 
Our  initial  estimate  from  figure  2a  for  the  first  even  solutions  was  n/2  -  d. 
For  each  successive  even  solution,  our  initial  guess  was  the  previous  final 
approximation  plus  n.  For  the  first  odd  solutions,  our  initial  guess  was 
k  -  d.  We  used  the  same  method  for  successive  solutions  (of  which  there 
were  none)  as  in  the  even  case.  The  markers  on  figure  2a  are  the  values  of 
the  solutions  generated  uy  the  program  with  the  Newton-Raphson  algo¬ 
rithm.  Their  "u"  value  is  given  in  table  1. 

We  now  wish  to  see  the  result  of  altering  our  input  parameters-thickness, 
dielectric  coefficient,  and  wavelength.  Table  2  shows  values  of  u  for  all 
possible  solutions  given  these  parameters.  All  combinations  were  used 
for  thicknesses  of  0.1,  0.35,  and  1.00  (am,  and  wavelengths  of  0.820, 
1.064,  and  1.550  Jim,  with  a  dielectric  coefficient  of  12.25.  We  also  did 
one  trial  with  a  dielectric  coefficient  of  4.80  (LiN03),  a  thickness  of  1 .00 
pm,  and  a  wavelength  of  1 .064  pm. 

We  can  see  from  this  table  that  increasing  the  wavelength,  decreasing  the 
thickness,  or  decreasing  the  dielectric  coefficient  will  decrease  the  number 
of  solutions.  However,  the  first  even  modes  (TE  and  TM)  will  never  dis¬ 
appear.  If  the  value  of  A  drops  below  Jt/2,  there  will  be  no  odd  solutions, 
and  only  one  of  each  even  solution.  If  the  value  of  A  is  very  high,  the 
number  of  solutions  will  be  approximately  4  x  A/k  ( A/n  of  each  type). 
The  separation  of  the  m’s  for  any  two  consecutive  solutions  of  the  same 
type  approaches  n  for  the  lower  value  modes.  These  results  are  illustrated 
in  figures  2a  through  j.  More  detailed  effects  of  altering  the  input 
parameters  are  shown  in  the  next  section. 
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igore  2.  Plots  of  both  sides  of  four  transcendental  equations  (19, 25, 33, 34)  presented  in  text.  Circles  mark 
points  derived  by  Newton-Raphson  approximation.  These  points  define  parameters  for  viable  modes  of 
propagation.  Calculations  were  done  for  a  slab  with  various  thicknesses,  dielectric  coefficients,  and  free 
space  wavelengths:  (a)  0.35  nm  slab  thickness,  12.25  dielectric  coefficient,  and  1.064-|im  free  space 
wavelength;  (b)  1.00-(im  slab  thickness,  4.80  dielectric  coefficient,  and  1.064-|im  free  space  wavelength; 

(c)  1.00-um  slab  thickness,  12.25  dielectric  coefficient,  and  1.550-(im  free  space  wavelength;  (d)  1.0-|im  slab 
thickness,  12.25  dielectric  coefficient,  and  1.064-fim  free  space  wavelength. 
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Figure  2.  Plots  of  both  sides  of  four  transcendental  equations  (J9, 25, 33, 34)  presented  in  text.  Circles  mark 
points  derived  by  Newton-Raphson  approximation.  These  points  define  parameters  for  viable  modes  of 
propagation.  Calculations  were  done  for  a  slab  with  various  thicknesses,  dielectric  coefficients,  and  free 
space  wavelengths  (cont’d):  (e)  1.00-pm  slab  thickness,  12.25  dielectric  coefficient,  and  0.820-pm  free  space 
wavelength;  (0  0.10-pm  slab  thickness,  12.25  dielectric  coefficient,  and  I.064-pm  free  space  wavelength; 

(g)  0.10-pm  slab  thickness,  12.25  dielectric  coefficient,  and  1.550-pm  free  space  wavelength;  (h)  0.10-pm  slab 
thickness,  12.25  dielectric  coefficient,  and  0.820-pm  free  space  wavelength. 


I? 


B 


f: 

vi 


ft 

Wl 


K 


HI 


j  Figure  2.  Plots  of  both  sides  of  Tour  transcendental  equations  (19, 25, 33, 34)  presented  in  text.  Circles  mark 

points  derived  by  Newton-Raphson  approximation.  These  points  define  parameters  for  viable  modes  of 
propagation.  Calculations  were  done  for  a  slab  with  various  thicknesses,  dielectric  coefficients,  and  free 
space  wavelengths  (cont’d):  (i)  0.35-)xm  slab  thickness,  12.25  dielectric  coefficient,  and  0.820>pun  free  space 
wavelength;  and  (j)  0  .35-pun  slab  thickness,  12.25  dielectric  coefficient,  and  1.550-pun  free  space  wavelength. 

i 

Table  1.  Values  of  u  for  all  possible  modes  with  slab  thickness  of 
0.35  pun  (d  =  0.175  |xm),  dielectric  constant  of  12.25,  and  free 
space  wavelength  of  1.064  pun 
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1.213211536 

3.373444796 


2.383396864 


1.53064(XX)6 

3.464896202 


3.001243353 
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4.  Calculation  of  Propagation  Parameters 


Once  we  have  found  values  for  u,  we  can  calculate  kv  T,  and  (3.  given  k 
13]  and  Ak  (the  bandwidth  of  the  source),  we  can  also  calculate  Ap  |4]  (the 
bandwidth  in  the  material).  We  shall  first  find  Ap  for  the  TE  even  mode. 
Differentiating  equations  (9),  (17),  (18),  and  (19)  gives  us 

k\  dk\  +  (3d(3  =  ekdk  (35) 

du  =dk\d  (36) 

dA  =  dkd<Tr\  (37) 


2  A(dA  du)  -  u  Va2  -  u2 

secz u  =  ,  -  = - 7 — 

tWA2  -u2  « 


Substituting  equation  (37)  for  dA  in  equation  (38)  and  rearranging  yield 


k  i  sec  u(u  tan  u  +  1) 


There  are  no  unknowns  on  the  right  side  of  the  equations.  We  can  now 
introduce  equations  (35)  and  (36)  and  solve  for  dp.  We  get 


kdk 

<*P  =  -j5~  Ve- 


e-  l 

sec2  u(u  tan  m  +  1) 


u  +  1)-^  ‘ 


For  the  other  solution  types,  we  get 


If  we  know  the  propagation  parameters  k,  T,  and  p,  we  can  find  expres¬ 
sions  for  £2,  £4,  H2,  and  H4  in  terms  of  Ev  £3,  Hv  and  Hy  Solving  equa¬ 
tion  (14)  for  £2  yields 
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£2  =£1  cos  k\derd  .  (41) 

Similar  manipulations  give  the  other  desired  relationships: 

£4  =  £3  sin  fci  er<i  (42) 

H2  =Hi  cos  Me™  (43) 

H4=Hisinkxderd  .  (44) 

These  equations  are  useful  for  calculating  the  electric  or  magnetic  field  at 
any  point,  which  will  be  done  later,  and  for  comparing  the  energy  prop¬ 
agated  inside  the  slab  to  that  outside  the  slab. 

The  time-averaged  energy  flux  at  any  point  is  defined  [5]  by 

<S>  =  /?e(^  Exh)  .  (45) 

In  the  TE  case,  Ez  is  defined  by  equation  (7)  inside  the  slab,  and  equation 
(8)  outside  the  slab.  Ex  and  Ey  are  both  zero.  Taking  H*  from  equations 
(3)  and  (4)  gives  us 


//*  =  ^  £1  cos  k\x  e 

\x\  <d 

(46) 

H*x  =  |  £2  e~rx  e-'fe 

M  >d 

(47) 

H)  =l-j-  Ex  sin k\xe~'*> 

\x\<d 

(48) 

H)  =  £  £2  e-rxe'^y 

1*1  >d  . 

(49) 

Inside  the  waveguide,  equation  (45)  becomes 

<S>  =  ^  £,  -|  cos2^|X  [x\<d  , 
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and  outside  we  have 
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We  define  the  efficiency  of  the  waveguide  as 

2  j  <S>  dx 

eff  =  ~z - 

2  f  <S>  dx 
o 

The  denominator  can  be  expressed  as 


2  j  <S>dx  =  2 J  <S>dx  +  2  j  <S>dx  . 

0  0  d 

We  first  calculate  the  integral  from  zero  to  d  (energy  inside  the  slab): 

2  /  <S>dx  =  f  E\  |  cos2  k\xdx 
o  471  K  o 


and  outside  the  slab  we  have 


2 1 <S>dx  =  ^£f22J e~2Vxdx  =  ^  <T2™  M>d.  (55) 


Now  using  equation  (41),  equation  (55)  becomes 

7  c|3£?  cos2 k\d 

2  j  <s>dx  =  Hrar- 

and  the  efficiency  of  the  waveguide  becomes 


e//  = 


Td  +  sin  u 

rd  +  1 
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for  the  TE  even  modes.  Similar  calculations  yield 


_ r "Y/  sin  u  /co\ 

eff~ — T3 — 

for  the  TE  odd  modes.  These  calculations  can  also  be  done  for  the  TM 
modes,  giving  us 


eff  = 


zVd  +  sin2  a 

TTd  +T' 


(59) 


for  TM  even  modes,  and 


*//  = 


er d  -  sin2« 

efd 


(60) 


for  the  TM  odd  modes. 

We  now  have  a  method  for  determining  the  propagation  parameters  of  a 
waveguide,  given  its  physical  properties  and  the  wavelength  used.  We 
can  also  determine  the  frequency  bandwidth  inside  the  material  for  each 
mode  of  operation,  given  the  bandwidth  of  the  incident  light.  The  effi¬ 
ciency  of  the  waveguide  is  also  calculable  for  each  mode.  These 
parameters  are  presented  in  table  3(a  through  j)  for  various  thicknesses, 
wavelengths,  and  indices  of  refraction. 

Also  presented  are  plots  of  Et  (or  Hz  for  the  TM  modes)  versus  x.  These 
plots  were  calculated  using,  for  the  TE  even  case,  equations  (7),  (8),  and 
(41),  and  assuming  that  £,  =  1.  For  the  other  cases,  we  used  correspond¬ 
ing  equations  and  assumed  either  £3  =  1,  Hx  =  1,  or  //3  =  1.  These  plots 
are  figures  3  through  8. 

An  interesting  feature  of  these  plots  is  the  bend  in  each  of  the  TM  curves. 
This  means  that  there  must  be  a  discontinuity  in  the  derivative  of  Hf.  This 
is  correct  and  is  caused  by  the  discontinous  polarization  current  density  at 
the  surface  of  the  slab.  In  calculations  for  the  TM  modes,  Ey  must  be  con¬ 
tinuous.  In  this  case,  Ey  is  defined  by 


At  x  =  d,  e  has  a  jump  discontinuity.  For  Ey  to  be  continuous,  d/dxH2 
must  also  have  a  jump  discontinuity. 

Interesting  conclusions  can  also  be  drawn  from  table  3.  We  see  that  A(J  is 
larger  than  Ak  and,  in  many  cases,  much  larger.  We  see  that  it  is  inversely 
proportional  to  p.  We  also  see  that  P,  as  expected,  can  have  a  value  from 
k  to  {e)mk.  We  can  also  compare  the  efficiencies  of  various  modes.  Sub¬ 
sequent  solutions  for  one  mode  type  are  always  less  efficient.  The  first 
TE  even  solution  is  always  more  efficient  than  the  second  TE  even  solu¬ 
tion,  and  so  on.  This  is  true  for  all  four  mode  types.  Also,  the  TM  modes 
are  always  more  efficient  than  their  corresponding  TE  modes. 

The  trends  evidenced  in  the  sample  data  cannot  necessarily  be  generalized 
for  all  situations. 

Table  3.  Propagation  Parameters  for  Slab  Waveguides 

Mode  ki  f  p  A$  efficiency 


A.  Thickness  =  1.00000  pm,  e  =  4.800,  wavelength  =  1.06400  pm, 
bandwidth  =  0.00100  pm,  k  =  59052,  and  A  k  =  55.501 


TEE 

26729 

111968 

126586 

123.47 

0.99183 

TEO 

53220 

102074 

117925 

138.91 

0.95812 

TME 

30281 

111060 

125784 

125.07 

0.99988 

TMO 

60284 

98067 

114474 

137.43 

0.99931 

TEE 

79097 

83636 

102382 

142.57 

0.90889 

TEO 

103363 

50669 

77811 

286.33 

0.68176 

TME 

89241 

72714 

93672 

167.92 

0.99667 

TMO 

111714 

27773 

65257 

241.68 

0.93811 

B.  Thickness  =  0.10000  pm,  e  =  12.250,  wavelength  =  0.82000  pm, 
bandwidth  =  0.00001  pm,  k  =  76624,  and  A  k  =  0.934 

TEE  169806  192919  207579  136  0.77780 

TME  249663  60993  97936  8.95  0.97879 


C.  Thickness  =  0.10000  pm,  e  =  12.250,  wavelength  =  1.06400  pm, 
bandwidth  =  0.00010  pm,  k  =  59052,  and  A  A  =  5.550 

TEE  146960  132792  145331  R23  0.66916 

TME  196598  24084  63775  62.91  0.87576 


Table  3.  Propagation  Parameters  for  Slab  Waveguide  (cont’d) 


Mode  ft, 

— r 

P 

AP 

efficiency 

D.  Thickness  =  0.10000  pm,  e 

=  12.250,  wavelength  =  1.55000  pm, 

bandwidth  =  0.00400  pm,  ft  =  < 

40537,  and  A  ft  = 

104.611 

TEE  1 14343 

73567 

83996 

324.79 

0.48295 

TME  135671 

8917 

41506 

1248.57 

0.60792 

E.  Thickness  =  0.35000  pm,  e  =  12.250,  wavelength  =  0.82000  pm, 
bandwidth  =  0.00001  pm,  ft  =  76624,  and  A  ft  =  0.934 

TEE 

73245 

246.347 

257989 

3.35 

0.98471 

TEO 

145209 

212052 

225472 

4.31 

0.91398 

88059 

241449 

253315 

3.46 

0.99998 

TMO 

175182 

188050 

203062 

4.32 

0.99986 

TEE 

143331 

162527 

4.42 

0.80362 

TME 

250183 

58822 

96599 

9.08 

0.99210 

F.  Thickness  =  0.35000  pm,  e  =  12.250,  wavelength  =  1.06400  pm, 
bandwidth  =  0.00010  pm,  ft  =  59052,  and  A  ft  =  5.550 


TEE 

69326 

185539 

194710 

20.07 

0.97115 

1 

U 

TEO 

136194 

143813 

155465 

33.22 

0.81213 

ft; 

TME 

87465 

177710 

187265 

21.44 

0.999% 

ft 

TMO 

171500 

99090 

115352 

34.81 

0.99908 

$ 

;«! 

TEE 

192768 

45512 

74556 

27.78 

0.47274 

ft 

i*’ 

TME 

197994 

5415 

59300 

67.53 

0.58392 

ft 

■ 

!•< 

J 

G.  Thickness  =  0.35000  pm,  e  = 

12.250,  wavelength  =  1.55000  pm, 

1 

1 

bandwidth  =  0.00400  pm,  ft  =  40537,  and  A  ft  = 

104.611 

i 

A 

TEE 

62477 

120759 

127382 

382.41 

0.93218 

TEO 

118791 

66143 

77577 

3650.06 

0.34053 

■ 

J 

TME 

85959 

105343 

112874 

460.23 

0.99981 

m 

TMO 

135486 

11391 

42107 

1236.25 

0.80131 

Table  3.  Propagation  Parameters  for  Slab  Waveguide  (cont’d) 


Mode  ^  T  P  Ap  efficiency 


H.  Thickness  =  1.00000  pm,  e  =  12.250,  wavelength  =  0.82000  (im, 
bandwidth  =  0.00001  tun,  *  =  76624,  and  A  *  =  0.934 


TEE 

29143 

TEO 

58258 

TME 

31216 

TMO 

62423 

TEE 

87315 

TEO 

1 16274 

TME 

93609 

TMO 

124758 

TEE 

145083 

TEO 

173657 

TME 

155836 

TMO 

186773 

TEE 

201845 

TEO 

229286 

TME 

217339 

TMO 

246044 

TEE 

254256 

TME 

256934 

255348 

266597 

250315 

261780 

255103 

266362 

249309 

260819 

241719 

253573 

229199 

241668 

239351 

251317 

224694 

2374(H) 

212139 

225553 

189460 

204368 

204369 

218262 

176544 

192455 

159092 

176583 

116102 

139108 

137170 

157120 

74256 

106701 

37494 

85306 

6041 

76862 

3.29 

0.99907 

3.36 

0.99589 

3.29 

l.(XXXX) 

3.36 

1.00000 

3.43 

0.99118 

3.69 

0.98214 

3.49 

0.99999 

3.69 

0.99999 

3.79 

0.97254 

4.50 

0.95180 

4.02 

0.99997 

4.56 

0.99993 

4.65 

0.93112 

7.26 

0.86289 

5.58 

0.99981 

8.22 

0.99850 

7.07 

0.65954 

11.40 

0.80353 

I.  Thickness  =  1.00000  (im,  e  =  12.250,  wavelength  =  1.06400  tun, 
bandwidth  =  0.00010  pm,  k  =  59052,  and  A  k  =  5.550 


TEE 

28526 

196003 

204706 

19.58 

0.99808 

TEO 

56994 

189691 

198670 

20.39 

0.99127 

TME 

31156 

195602 

204322 

19.65 

1.00000 

TMO 

62291 

188018 

197074 

20.37 

0.99999 

TEE 

85339 

178741 

188243 

20.96 

0.98132 

TEO 

113464 

162348 

172754 

24.22 

0.95957 

TME 

93376 

174677 

184389 

21.77 

0.99998 

TMO 

124349 

154170 

165093 

24.32 

0.99995 

TEE 

141207 

138894 

150927 

25.04 

0.93603 

TEO 

168205 

104585 

120105 

38.66 

0.86209 

TME 

155033 

123270 

136685 

29.37 

0.99986 

TMO 

184392 

72323 

93370 

43.00 

0.99906 

TEE 

193020 

44431 

73900 

39.62 

0.70521 

TME 

197898 

8205 

59620 

67.29 

0.86810 

Table  3.  Propagation  Parameters  for  Slab  Waveguide  (cont’d) 


Mode 

*« 

— r — 

P 

AP 

efficiency 

J.  Thickness  =  1.00000  pm,  e  =  12.250,  wavelength  =  1.55000  pm, 
bandwidth  =  0.00400  pm,  k  =  40537,  and  A  A c=  104.61 1 

TEE 

27363 

133182 

139215 

371.33 

0.99471 

TEO 

54572 

124532 

130963 

407.88 

0.97413 

TME 

31033 

132375 

138443 

375.23 

l.(XXKX) 

TMO 

61996 

121007 

127617 

407.06 

0.99998 

TEE 

81409 

108898 

116198 

424.22 

0.94437 

TEO 

107441 

83323 

92661 

662.16 

0.85012 

TME 

92728 

99437 

107382 

483.76 

0.99991 

TMO 

122329 

59344 

71868 

722.82 

0.99924 

TEE 

131055 

36207 

54352 

665.58 

0.66940 

TME 

135824 

6171 

41004 

1265.66 

0.84028 

Figure  3.  Electric  field 
(arbitrary  units)  versas 
distance  from  slab  center 
for  two  even  and  two  odd 
TE  modes,  for  a  slab  of 
1  pm  thickness,  a 
dielectric  coefficient  of 
4.80,  and  a  free  space 
wavelength  of  1.064  pm. 
Vertical  line  in  this,  and 
all  following  figures,  rep¬ 
resents  slab  edge. 
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Figure  6.  Klectric  field 

(c) 

o 

versus  distance  from  slab 

center.  Slab  thickness  of 

035  pm,  a  dielectric 

coefficient  of  12.25,  and 

a  m 

free  space  wavelength  of 

>N 

<0 

(a)  1.550  pm,  for  one  even 

n 

mode  and  one  odd  mode; 

1  9 
<  o 

(b)  1.064  pm,  for  two  even 

"O 

modes  and  one  odd  mode; 

.2 

and  (c)  0.820  pm,  for  two 

is 

even  modes  and  one  odd 

%  « 

mode. 

UJ 
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Figure  7.  Magnetic  field 
versus  distance  from  slab 
center.  Slab  has  thickness 
of  0.35  pm,  a  dielectric 
coefficient  of  12.25,  and 
free  space  wavelength  of 

(a)  1.550  pm,  for  one  even 
mode  and  one  odd  mode; 

(b)  1.064  pm,  for  two  even 
modes  and  one  odd  mode; 
and  (c)  03120  pm,  for  two 
even  modes  and  one  odd 
mode. 
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C$FORT  SLAB1 
CfLINK  SLAB1 
CfPURGE  SLAB1 . EXE 
CfDEL  SLABl . OBJ ; • 

C 

C  THIS  PROGRAM  FINDS  THE  SOLUTIONS  FOR  AN 

C  INFINITE  PLANAR  SLAB  WAVEGUIDE  WITH  DIELECTRIC 

C  CONSTANT  EPS,  THICKNESS  THK ,  AND  INCIDENT 

C  WAVE  LENGTH  WL. 

C 

PROGRAM  NEWRAP 
CHARACTER* 1  ANS 

DIMENSION  TEE ( 500 ) ,TEO(500) ,TME(500) ,TM0(50C) 

TYPE  *,•  THIS  PROGRAM  SOLVES  FOR  THE' 

TYPE*,'  ROOTC  TO  0'rP.  SLAB  WAVEGUIDE  PROBLEM.' 

TYPE* , '  ' 

C  TYPE*,'  ENTER  DIELECTRIC  CONSTANT , TH I CKNESS , WAVE  LENGTH' 

C 

C  THIS  FILE  CONTAINS  THE  RUN  PARAMETERS.  IT  IS  MADE  BY  THE  USER. 

C  YOU  MUST  STORE  THEM  IN  THE  FILE  LIKE  THIS: 

C  EPSILON  THICKNESS (cm)  WAVELENGTH (cm)  BANDWIDTH) cm) 

C 

OPEN ( UN I T=9 , TYPE= ’ OLD ' . NAME= ' SLABRUN . DAT ' ) 

READ ( 9 , * )  ,  EPS  ,  THK.WL.BW 
THK=THK/2 
C 

C  SET  RUN  CONSTANTS 

C 

EFS=l .0 

PI=3. 141592554 

A=2*PI *THK* (EPS-1 ) ** . 5/WL 

IO=IFIX( A/PI+ . 5 ) 

C 

C  IF  A  IS  LESS  THEN  PI/2  THEN  WE  HAVE  ONLY  TWO  ROOTS,  ONE  TE  EVEN 

C  AND  ONE  TM  EVEN.  THERE  WILL  BE  NO  ODD  ROOTS,  SO  WE  GO  TO  THE 

C  SPECIAL  EVEN  ONLY  ROUTINE. 

C 

IF( A . LT . PI/2)GOTO  123 
100  FORMAT) 16X.4F16. 9) 

C 

C  START  LOOP.  AND  SET  INITIAL  GUESSES  FOR  THE  NEWTON-RAPHSON 

C  ROUTINE 

C 

1  =  1 

SE=P I /2 - O-THK 
SE1=SE 
SO=P I -THK 
C 

C  IF  OUR  GUESS  IS  GREATER  THAN  A,  THE  PROGRAM  WILL  CRASH.  SO  WE 

C  INSURE  THAT  OUR  GUESSES  ARE  LESS  THAN  A. 

IF  (SO  .GT.  A ) SO=A-THK 
S01=S0 
C 

C  USE  SUBROUTINES  TO  FIND  ROOTS 

C  EPS  IS  EPSILON,  EFS  IS  EPSILON  OF  FREE  SPACE 

C 

C  THERE  SHOULD  BE  ONE  SOLUTION  OF  EACH  TYPE  -  TWO  ODD,  TWO  EVEN, 

C  OF  TM  AND  TE  MODES. 

C 

10  IF  (SO  .GE.  A)  SO=A-THK 

IF  (SOI  .GE.  A)  S01=S0 
CALL  ODD(TEO( I ) , SO, EFS, A) 

CALL  ODD) TMO( I ) , SOI , EPS , A ) 

CALL  EVEN(TEE( I ) ,SE,EFS, A) 

CALL  EVEN(TME( I ) ,SE1 ,EPS,A) 

C 

C  SET  UPPER  AND  LOWER  LIMITS  FOR  EACH  INDIVIDUAL  MODE 
C 

AI=PI*(2*I+l)/2.0 
AJ=PI *( 2*1-1 )/2 . 0 
AL=Pl*(2*I-3)/2.0 
C 


CHECK  THAT  ROOT  IS  WITHIN  THE  SET  LIMITS 

THAT  IS,  CHECK  THAT  NO  ROOTS  ARE  SKIPPED,  OR  REPEATED. 

IF  (TEO(I)  .GT.  A I  .OR.  TMO( I )  .GT.  AI)  1 FLG= I FLG+ 1 

IF  (TEE(I)  .GT.  AJ  OR.  TME(I)  .GT.  A J)  JFLG=JFLG+1 

IF  (TEO(I)  .LT.  AJ  OR.  TMO( I )  .LT.  AJ )  KFLG=KFLG+1 

IF  (TEE(I)  .LT.  AL  OR.  TME( I )  .LT.  AL)  LFLG=LFLG+1 

TYPE  100 ,TEE( I ) ,TEO( I ) ,TME( I ) ,TMO( I ) 

1  =  1  +  1 

CHECK  TO  SEE  IF  WE'RE  DONE. 

IF  (I  .LE.  10)  GOTO  10 

CHECK  TO  SEE  IF  WE  HAVE  ANOTHER  PAIR  OF  EVEN  ROOTS. 

IF  ( A— ( I — 1 ) • P I  .LE.  0)  GOTO  987 

FIND  PAIR  OF  EVEN  ROOTS 

SE=A-THK 

SEl=SE+THK/2 

TYPE* . I0*PI , SE, A 

CALL  EVEN(TEE( I ) ,SE,EFS,A) 

CALL  EVEN ( TME ( 1 ) , SE1 , EPS , A ) 

IXFL=1 

IF  ANY  OF  THE  FLAGS  BELOW  ARE  NOT  ZERO,  A  MODE  WAS  SKIPPED 
OR  REPEATED, 

TYPE  * , I FLG , JFLG , KFLG , LFLG 
FORMAT ( A 1 ) 

OPEN  FILE  FOR  EZ-GRAPH  PLOTTING,  AND  FILE  FOR  PRINT  OUT 

OPEN ( UN I T=2 . TYPE= ' NEW ' , NAME= ’ ALROOT . EZG ’ ) 

OPEN  (UNIT=1 , TYPE= ’ NEW ' , NAME= ' ALROOT . DAT ' ) 

WRITE(1,*)  ’  ' 

DO  20  1=1  ,  10,1 

WRITE( 1 . 100)  TEE( I ) ,TEO( I ) ,TME( I ) ,TMO< I ) 

CALCULATE  THE  VALUES  FOR  THE  PLOTTING  OF  THE  INTERSECTION 
OF  THE  TRANSCENDENTAL  EQUATIONS. 

CN=-1 . O/TAN ( TEO( I ) ) 

CN 1 =- 1 . 0/TAN(TM0( I ) ) 

TN=TAN(TEE( I ) ) 

TN1=TAN(TME( I ) ) 

WRITE  TO  EZG  FILE 

WRITE( 2 , * )  TEE( I ) ,TN 
WRITE( 2 , * )  TEO( I ) ,CN 
WRITE( 2 , * )  TME( I ) ,TN1 
WRITE(2 , * )  TMO( I ) ,CN1 
CONTINUE 

IF  ( IXFL  NE.  1)  GOTO  789 
WR!TE(2, • )  0.0,99999. 

WRITE( 1,234)  TEE( I ) , TME( I ) 

FORMAT(  16X  ,F16 . 9,  '  '.Fie.D) 

TN=TAN(TEE( I ) ) 

WRITE( 2 , • )  TEE( I ) ,TN 
TN=TAN(TME( I ) ) 

WRITE( 2, • )  0.0,99999. 

WRITE( 2 , * )  TME( I ) ,TN 

STOP 

END 

SUBROUTINE  ODD(F,X,E,A) 

THIS  SUBROUTINE  FINDS  THE  ROOTS  OF  ODD  SOLUTION  MODES 
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DEPENDING  ON  E,  IT  WILL  SOLVE  FOR  TWO,  OR  TME 

DOUBLE  PRECISION  Q.F2.F3.D1 
Q=l/(A*A-X*X)**.5 
TYPE* ,X,D1 
F2=TAN(X)+X*<J/E 

F3=C0S ( X ) *  » ( -2 ) +Q/E+X  *  X • Q  *  *  3/E 
Dl=-F2/F3 

IF  (ABS(Dl)  .LE.  .00001 )GOTO  310 

X=X+D1 

GOTO  320 

F=X 

X=X+3. 141692554 

RETURN 

END 

SUBROUTINE  EVEN ( F , X , E , A ) 

THIS  SUR0UTINE  FINDS  THE  ROOTS  FOR  EVEN  SOLUTION  MODES 

DOUBLE  PRECISION  Q,F2,F3,D1 

Q=(A*A-X*X)**.5 

TYPE* ,X,D1 

F2=TAN ( X ) -Q*  E/X 

F3=COS(X) *  * (-2 )+Q*E/(X*X)+E/Q 

Dl=-F2/F3 

IF  (ABS(Dl)  LE.  .00001)  GOTO  410 

X=D1+X 

GOTO  420 

F=X 

X=X+3. 141592554 

RETURN 

END 
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C*FORT  SLAB2 
CtLINK  SLAB2 
C$DEL  SLAB2 . OBJ ; • 

CtPURGE  SLAB 2 . EXE 
C 

C  THIS  PROGRAM  CALCULATES  THE  PROPAGATION  PARAMETERS  FOR  THE 

C  INDIVIDUAL  SOLUTIONS  TO  THE  SLAB  WAVEGUIDE  PROBLEM. 

PROGRAM  BETAWD 
CHARACTER* 5  MD(4) 

REAL  K ,Kl ,EFF(4) 

DATA  PI/3.1415926538/ 

C 

C  THESE  CHARACTER  STRINGS  REPRESENT  THE  FOUR  DIFFERENT  MODE  TYPES. 
C 

MD( 1  )  = '  TEE  1 
MD(2)='  TEO 
MD(3)=’  TME  1 
MD( 4 )= '  TMO  ■ 

C 

C  OPEN  THE  RUN  PARAMETER  FILE (CREATED  BY  THE  USER). 

OPEN ( UN  I T=9 . TYPE= ' OLD ’ , NAME= ' SLABRUN . DAT ' ) 

C 

READ( 9 . * )  EPS , THK . WL , DW 
DK=2»PI*DW/(WL**2-DW**2) 

D=THK/2 

K=2*PI/WL 

C 

C  OPEN  SOLUTION  DATA  FILE,  CREATED  BY  SLAB1  PROGRAM) 

OPEN  ( UN I T=2 , TYPE= ' OLD ' , NAME=  * ALROOT . EZG ' ) 

C 

C  OPEN  OUTPUT  FILE  FOR  PROPAGATION  PARAMETERS,  AND 

C  WRITE  FILE  HEADER 

OPEN  ( UN I T= 1 , TYPE= ' NEW ' , NAME= 1 PARAMS . DAT • ) 

WRITE(1,*)'  Table  .  Propogation  parameters  for  slab  waveguide 
+  with' 


WRITE(1,*) 

WRITE(l.lO) 

THK *10000, EPS 

10 

FORMAT ( ' 

thickness=' ,F7.5,  '  urn' 

,  ’  eps i lon=- ,F7 . 

WRITE( 1,12) 

WL* 1 0000, DW* 10000 

12 

FORMAT ( ' 

+F7 . 5 , '  urn ’ ) 

wavelength=' ,F7.5, '  um' 

, ‘  bandwidth=' , 

WRITE( 1,14) 

K.dk 

14 

FORMAT ( ' 
wr  i  te( 1 , * ) ' 

k= ' , F7 . 0, ' 

delta  k=' ,F7.3) 

+ 

• 

WRITE(1,*)' 
WRITE(1 ,*) 

'mode  k  gamma 

beta 

+  delta  beta 

ef f iciency 

wr  ite( 1 ,*) 
writed  ,*) ' 

•  1  ’ 

+. 


999  FORMAT (A5 ,3F10.0,F11 .2.F13.5) 

C 

C  MODE  TYPE  IS  KEPT  TRACK  OF  WITH  MTP(0  BEFORE  PROGRAM  STARTS) 

MTP=0 

100  READ( 2 , • , END=777 )  U, DUMMY 

MTP=MTP+ 1 

IF  (MTP  .GT.  4)  MTP=MTP-4 
I F ( U . EQ .  0)  GOTO  100 
C 

C  CALCULATE  K1 .GAMMA 

K1=U/D 

GAM=(K*K*( EPS-1 )-Kl*Kl)**.5 
C 

C  CALCULATE  EFFICIENCY  AND  BETA  "B" 

B=(GAM**2+K*K)**.5 

SU=SIN(U) 

CU=COS(U) 

sc=su*cu 

GD=GAM*D 

GED=GD*EPS 

EFF( 1 )=(GD+SU*  *2 )/( GD+1 ) 


EFF ( 2 ) = ( GD-SU • • 2 ) /GD 
EFF(  3  )=(GED+SU**2 )/(GED+l  ) 

EFF ( 4 ) = ( GED-SU • • 2 ) /GED 

CALCULATE  THE  BANDWIDTH  OF  BETA 

CSU2=l/SU**2 

SCU2=l/CU**2 

tu=su/cu 

ct=l/tu 

Q=1 

IF  (MTP  .GT.  2)  Q=l/EPS*  *2 

if  (mtp  .eq.  2* int (mtp/2 ) )  goto  987 

t 1=9cu2 

f  2=tu 

goto  986 

f 1=c9u2 

f  2=et 

db=k*dk*(eps-q*(ep9-l )/(fl*(u*f2+l)))/b 
WHITE  TO  OUTPUT  FILE 

WKITE(1 ,999)  MD(MTP) ,K1 , GAM , B , DB , EFF ( MTP ) 
GOTO  100 

wr  i  te( 1 ,*)■ _ 


C#FORT  SLAB3 
CfLINK  SLAB 3 
C$DEL  SLAB3 . OBJ ; • 

CfPURGE  SLAB3.EXE 
C 

C  THIS  PROGRAM  CALCULATES  THE  VALUES  OF  THE  CURVES  REPRESENTED 

C  BY  THE  FOUR  TRANSCENDENTAL  EQUATIONS  FOR  OUR  SLAB  WAVEGUIDE 

C  PROBLEM.  THE  FOUR  CURVES  NEEDED  ARE. 

C 

C  FI (U)=TAN(U)  F2(U)=-C0T(U)  F3(U )  =  ( A»A-U*U) *  * . 5  F4(U)=EPS»F(3) 

C 

PROGRAM  RTGRA 
DIMENSION  F(  4 ) 

DATA  U ,  WL ,PI/0 . 0 , .000106,3.1415926536/ 

C 

C  OPEN  RUN  PARAMETFR  FILE  (CREATED  BY  THE  USER)  AND  READ  PARAMETERS 
OPEN ( UN  I T=9 , TYPE= ' OLD ' , NAME= 1 SLABRUN . DAT ' ) 

READ( 9 , • )  EPS , THK , WL , BW 
C 

C  CALCULATE  D(HALF  THICKNESS)  AND  A (DESCRIBED  IN  PAPER) 

D=THK/2 

A=2*PI *D* (EPS-1 ) • • . 5/WL 
N=0 

DU-P 1/2000. 

C 

C  OPEN  FILE  FOR  OUTPUT  OF  CURVES.  THE  OUTPUT  IS  OF  THE  FORM: 

C  U  F( 1 )  F ( 2 )  F(3)  F( 4 ) 

C 

9  0PEN(UNIT=1 ,TYPE='NEW' , NAME= ' SMROOT . DAT ' ) 

C 

C  BEGIN  CURVE  CALCULATING  LOOP 

DO  10  1=0,999 
C 

C  WE  CALCULATE  IN  BLOCKS  OF  PI/2  TO  AVOID  A  VERTICLE  LINE  WHEN 

C  WE  PLOT  THE  TRIGONOMETRIC  CURVES  WITH  EZGRAPH. 

U=I *DU+N*PI/2 . 

IF(U  .GT.  A)  GOTO  11 

IF  ( INT( 2*U/PI )  .EQ.  2*U/PI )  GOTO  10 

F( 1 )=TAN(U) 

F(2)=-l/F( 1 ) 

F(3)=( A*A-U*U) •• . 5/U 
F(4)=EPS*F(3) 

DO  20  J=1  ,4 

IF( ABS(F( J ) )  .LT.  100)  GOTO  20 
C 

C  FOR  PLOTTING  PURPOSES,  WE  DONT  WANT  ANY  OF  THE  FUNCTIONS  VERY  HIGH 
C  IT  DOESN'T  MATTER  MUCH  EITHER  WAY 
F( J )=99 . 9999*ABS( F( J ) )/F( J  ) 

20  CONTINUE 

WRITE( 1,100)  U,F(l),F(2),F(3),F(4) 

10  CONTINUE 

100  FORMAT( 5F15 . 10) 

N=N+1 

11  CL0SE(UNIT=1) 

IF(U  LT.  A  )  GOTO  9 
STOP 


C$FORT  SLAB4 
C#L INK  SLAB4 
CfDEL  SLAB4 . OBJ ; * 

ClPURGE  SLAB4.EXE 
C 

C  THIS  PROGRAM  CALCULATES  THE  VALUES  OF  THE  ELECTRIC 

C  OR  MAGNETIC  FIELDS  FOR  ANY  X,  INSIDE  OR  OUTSIDE  THE 

C  SLAB.  IT  ASSUMES  A  MAXIMUM  POSSIBLE  FIELD  OF  ONE,  OR 

C  El=l ,  AND  Y=2*PI *n . 

PROGRAM  MODEPL 
REAL  K,K1 

DATA  PI/3.1415926536/ 

COMMON  U.D.K.Kl , B , GAM , EGD , EPS ,CU , SU , EO 
C 

C  OPEN  FILE  CREATED  BY  USER  WITH  RUN  PARAMETERS 
OPEN ( UN  I T=9 , TYPE= ' OLD ’ ,NAME= ’ SLABRUN .DAT' ) 

READ( 9 , • )  EPS ,THK ,WL ,BW 
C 

C  CALCULATE  PARAMETERS  USED  FOR  FINDING  FIELDS 
D=THK/2 
K=2*P1/WL 
C 

C  OPEN  FILE  WITH  SOLUTIONS  TO  TRANSCENDENTALS .CONTAIN I NG  VALUES  FOR  U 

C  THIS  FILE  WAS  CREATED  BY  SLAB1 

OPEN  ( UNI T=2 , TYPE= ' OLD ' , NAME= ' ALROOT . EZG ' ) 

C 

C  KEEP  TRACK  OF  MODE  TYPE  WITH  MTP 
MTP=-1 

100  READ ( 2 , • , END=777 )  U, DUMMY 

MTP=MTP+ 1 

IF(U.EQ.  0)  GOTO  100 
CU=COS(U) 

SU=SIN(U) 

C 

C  CALCULATE  PROPAGATION  PARAMETERS 

K1=U/D 

GAM=( K*K* ( EPS-1 )-Kl*Kl)**.5 
B=(GAM**2+K*K)**  .  5 
EGD=EXP(GAM*D) 

C 

C  OPEN  FILE  FOR  OUTPUT 

0PEN(UNIT=1 , TYPE= ' NEW ’ , NAME= ' SMMODE . DAT ' ) 

DX=10 . "D/1000 . 

C 

C  LOOP  TO  CALCULATE  FIELDS 
DO  200  1=0,1000 
X=I *DX 

E0=MTP-2* I F I X ( FLOAT ( MTP ) /2 . 0 ) 

CALL  MODCLC ( X , V ) 

WRITE( 1 , • )X*10000,V 
200  CONTINUE 

CL0SE(UNIT=1 ) 

GOTO  100 
777  STOP 

END 

SUBROUTINE  MODCLC (X,V) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  FIELD  DEPENDING  ON  WHICH  TYPE  OF 
C  IS  BEING  PLOTTED 
REAL  K.Kl 

COMMON  U,D,K,K1 ,B .GAM , EGD , EPS ,CU , SU ,E0 
C 

C  CHECK  FOR  INSIDE  OR  OUTSIDE  OF  SLAB 
IF  (X  .LT.  D)  GOTO  9 


CHECK  FOR  EVEN  OR  ODD  MODE 
IF( EO)  2,2,1 
AM=CU 
GOTO  3 
AM=SU 

V=AM»EGD»EXP(-X*GAM) 
GOTO  4 

IF  (EO)  6,8,6 
V=COS(Kl»X) 

GOTO  4 
V=SIN(K1*X) 

RETURN 
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